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Abstract 

A model is derived for the viscoelastic response of glasses at isothermal uniaxial 
deformation with small strains. A glass is treated as an ensemble of relaxing units 
with various activation energies for rearrangement. With reference to the energy- 
landscape concept, the rearrangement process is thought of as a series of hops of 
relaxing units (trapped in their potential wells on the energy landscape) to higher 
energy levels. Stress-strain relations are developed by using the laws of thermody- 
namics. Adjustable parameters are found by fitting experimental data in torsional 
dynamic tests on a multicomponent silicate glass at several temperatures near the 
glass transition point. 
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1 Introduction 

This study is concerned with the viscoelastic behavior of silicate glasses at uniaxial defor- 
mation with small strains. This subject has been a focus of attention in the past decade 
because of numerous applications of vitreous silica in industry (that range from waveg- 
uides and optical fibers to semiconductor wafers and tire additives |J and the key role of 
multicomponent silicate glasses in geological processes [0. 

Observations in static and dynamic mechanical tests on silicate glasses are conven- 
tionally fitted by using such phenomenological approaches as (i) the generalized Maxwell 
model ||, (ii) the KWW (Kolhrausch-Williams-Watt) formula ||, and (iii) linear stress- 
strain relations with fractional derivatives ||. The objective of this paper is to develop 
constitutive equations for the linear viscoelastic response of an amorphous glass based 
on a micro-mechanical concept and to validate these relations by matching experimental 
data in torsional oscillatory tests on a multicomponent silicate glass. 
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2 A micro- mechanical model 



A silicate glass is thought of as an ensemble of cooperatively rearranging regions (CRRs), 
where mechanical stress relaxes due to rotational rearrangement of small clusters of atoms 
(consisting of a few SiC>4 tetrahedra) that change their configuration as they are agitated 
by thermal fluctuations. 

To describe the viscoelastic response of a glass in dynamic tests with a restricted win- 
dow of frequencies, we distinguish two types of relaxing units: (i) active CRRs that can 
rearrange during the experimental time-scale, and (ii) passive CRRs whose characteristic 
time for rearrangement substantially exceeds the duration of a test. An active CRR is en- 
tirely characterized by the energy, v, of thermal fluctuations necessary for rearrangement. 
With reference to the energy-landscape concept , an active CRR is modelled as a mate- 
rial point trapped in a potential well on the energy hypersurface in the phase space. The 
well is characterized by its depth, v, with respect to some reference (liquid-like) energy 
level. The point spends most time being located at the bottom level of its potential well. 
At random instants, the point hops to higher energy levels as it is thermally agitated. If 
the point does not reach the liquid-like energy level in a hop, it returns to the bottom 
level of its well being unaltered, when the point reaches the reference level in a hop, the 
CRR changes its configuration and accepts the current state of a deformed medium as its 
reference state. 

The kinetics of rearrangement is uniquely determined by the temperature-dependent 
attempt rate To (the number of hops in a well per unit time) and the distribution function, 
q(u), for hops with various heights (q(u)du is the probability of a hop whose height belongs 
to the interval [u, u + du]). We adopt the exponential distribution of intensities of hops 

?(") = y eX P( _ y) ( M >°)> 

where V is the average height to be reached in a hop. For a CRR trapped in a well with 
depth v, the probability of rearrangement in a hop reads 

roo v 

Q( v ) = J q(u)du = exp(-— ). 

The rate of rearrangement, 

T(v)=T Q(v), 

equals the number of hops (per unit time) in which a point reaches the liquid-like level, 

r» = r ex P (-^). (i) 

Denote by iV a and N p the numbers of active and passive CRRs per unit mass, and by p(v) 
the distribution of potential wells with various energies v. The rearrangement process is 
described by the function n(t, r, v) that equals the number (per unit mass) of active CRRs 
trapped in wells with energy v which have last been rearranged before instant r e [0, t]. 
The quantity n(t,t,v) is the current number of active CRRs (per unit mass) with energy 
v, 

n(t,t,v) = N a p(v). (2) 
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The quantity 



dn 

— {t,T,V)dT 



equals the number (per unit mass) of active CRRs with energy v that have last been 
rearranged during the interval [r, r + rfr] and, afterwards, have not beed rearranged within 
the interval [r, t] . The number of active CRRs (per unit mass) with energy v that were 
rearranged (for the first time) within the interval [t, t + dt] is given by 

-—(t,0,v)dt, 

and the number (per unit mass) of active CRRs with energy v that have been rearranged 
within the interval [r, r+dr], and, afterwards, were rearranged during the interval [t, t+dt] 
reads 

O^n 

The rate of rearrangement, T, is defined as the ratio of the number of active CRRs 
rearranged per unit time to the total number of active CRRs, 

on 

^(t,0,v) = -l»ra(t,0,u), 

d 2 n dn 
g^(t,r,v) = -r(v)-(t,r,v). 

The solutions of these equations with initial condition (2) read 



n(t,0,v) = N a p(v) exp -F(v)t 

— (t,r,v) = N aP (v)r(v)ew[-r(v)(t-T) 



(3) 



3 Constitutive equations 

At uniaxial deformation with small strains, a CRR is modelled as a linear elastic medium 
with the strain energy 

1 2 

w = -lie , 

where \i is an average rigidity per CRR and e is the strain from the reference state of a 
relaxing unit to its deformed state. We assume that CRRs are connected by links that 
transmit the macro-strain in an ensemble, e, to individual regions. This means that for 
a passive CRR (that is not rearranged within the experimental time-scale), as well as for 
an active CRR that has not been rearranged before time t, the strain e coincides with 
the macro-strain e(t). For a CRR that has last been rearranged at time r e [0,t], the 
reference state coincides with the state of the ensemble at the instant of rearrangement, 
which implies that the strain, e, equals the difference between the current macro-strain, 
e(t), and the macro-strain, e(r), at the instant of rearrangement, r. 

To calculate the strain energy density per unit mass, W(t), we multiply the mechanical 
energy per CRR by the number of relaxing units and sum the mechanical energies of 
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passive CRRs and active CRRs rearranged at various instants r e [0,t\. Neglecting the 
energy of interaction between relaxing units, we find that 



W{t) = -p 



N p + n(t,0,v)dv 

J 



e\t) 



t dn 



+ J o dv j —(t,T,v)[e(t)-e(r) 

Differentiation of Eq. (4) with respect to time results in 

dW , s . . s de , s , r , s 
— { t ) = A( t )-(t)-Y(t), 



'dr\. 



A(t) = pN 
1 



:(*) - k j\{r)dT J™ T{v) exp(-r»(£ - rj)p(v)di 



Y(t) 



-J' 



F(v)dv 



rt , \ 2 

n{t,0,v)e 2 {t) + J —(t,r,v)(e(t)-e(r)) dr 



(4) 

(5) 
(6) 
(7) 



where N = N a + N p and k = N a /N. For isothermal uniaxial deformation, the Clausius- 
Duhem inequality reads 

dW a de ^ ^ 
dt p dt ~ ' 

where p is mass density, a is stress, and 9 is internal dissipation per unit mass. Substi- 
tution of Eq. (5) into this equation yields 

= -(a- P A)^ + Y > 0. 
p dt 

It follows from Eq. (7) that Y(t) is non-negative. This implies that the dissipation in- 
equality is satisfied, provided that the first term on the right-hand side vanishes. Equating 
the expression in brackets to zero and using Eq. (6), we arrive at the constitutive equation 



a(t) = G 



s(t) - k J* e{r)dT J™ F(v) exp(-F(v){t - r))p{v)di 



(8) 



where G = ppN is an elastic modulus. 



4 Dynamic mechanical tests 

For an oscillatory test with an amplitude e and a frequency u, the strain reads 

e(t) = €o exp(iut), 

where % = T- Substituting this expression into Eq. (8) and using Eq. (1), we find the 
storage, G'(lu), and loss, G"{u), moduli 

G ' M = s ['-< ii^ ^]. 
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where 

v 

Z ~yi K = kV i P{ z ) = P(Vz). 

Restrictions on the upper boundary of a frequency window in conventional mechanical 
tests (the maximal frequency does not exceed 1 kHz) imply that the shape of the distri- 
bution function P{z) at small z cannot be uniquely determined by fitting observations. 
To avoid ambiguities, we set 

P(z) = K 1 5(z) + K 2 P (z), (10) 

where 8(z) is the Dirac delta-function, K m (m = 1,2) are non-negative constants such 
that Ki + K 2 = 1, and Pq{z) is the distribution function corresponding to the range 
of frequencies "visible" in the frequency window. With reference to ||, the following 
expression is adopted for the function Pq{z): 

Po{z) = 6z(z) 6XP ("z)' 

where 
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and (z) is the average activation energy per CRR. Substitution of Eq. (10) into Eq. (9) 
results in 

g\u) = g u " 2 ■ - ■ - • - 2 r Poiz)dz 



G"(u) = Geo 



T^ + uj 2 Jo T 2 exp(-2z) +u 

T f<x>r Q exp(-z)P Q {z)dz 



Tl + u 2 Jo T 2 exp(-2z) +u 2 



(12) 



where ki = KqKi, k 2 = K0K2, and k 3 = 1 — Kq. The coefficient k\ is the concentration 
of CRRs whose rates of rearrangement substantially exceeds the maximal frequency of 
oscillations ("viscous" units), k 2 is the concentration of CRRs whose rates of rearrange- 
ment belong to the frequency window of dynamic tests ( "viscoelastic" units), and k% is 
the fraction of CRRs whose rates of rearrangement is noticeably less than the minimal 
frequency of vibrations ("elastic" units). 

Equations (11) and (12) are determined by 5 material constants: 

1. the elastic modulus G, 

2. the attempt rate r , 

3. the average activation energy per CRR (z), 

4. the fraction of "viscous" CRRs k\, 

5. the fraction of "elastic" CRRs k%. 

These parameters are found by fitting experimental data for a multicomponent silicate 
glass. 
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5 Validation of the model 



We study the mechanical response of a rock wool glass with the chemical composition (wt.- 
%) 41.9 Si0 2 , 17.0 A1 2 3 , 18.3 CaO, 12.5 MgO, 0.3 K 2 0, 1.4 Ti0 2 , 5.7 FeO, 2.0 Na 2 0. 
Torsional dynamic tests were performed by N.S. Bagdasarov (University of Frankfurt, 
Germany) at temperatures ranged from 500 to 760 °C (this interval includes the glass 
transition point T g = 678 °C). For a description of specimens and the experimental 
procedure, see 0. The effect of temperature on the viscielastic behavior of the glass is 
demonstrated by Figure 1, where the loss tangent (at the frequency / = 1Hz) is plotted 
versus temperature. To approximate observations, we present Eqs. (12) in the form 



CO 



Ci 



Po(z)dz 



c 1 r 



+ c 2 



r^exp(-2^) 
r exp(—z)Po(z)dz 
rgexp(-2^) +cu 2 



(13) 



where C„ 



Gk r , 



1,2,3). 



To find r , Z and C m at a given temperature, T, we fix some intervals [0, r max ] 
and [0, Z max ], where the "best-fit" parameters r and Z are assumed to be located, and 
divide these intervals into J subintervals by the points FV> = jAT and Z^ = kAZ 
(J, k = 1, . . . , J) with Ar = T max /J and AZ = Z mSLX /J. For any pair, {T^\ Z<*)}, the 
constants C m = C m (j, k) are found by the least-squares technique from the condition of 
minimum of the function 



M(j,k) 



E 



[( G cxp( 



— G r 



where G' ew and G'^ p are the storage and loss moduli measured in a test, the quanti- 
ties G' nnm and G" um are given by Eqs. (13), and the sum is calculated over all fre- 
quencies The "best-fit" parameters To and Z minimize the function M on the set 
{rO), Z^ (j,k = 1, . . . , J)}. After the constants C m are determined, the elastic modu- 
lus, G, and the constants k m are calculated by the formulas G = Ci+C 2 +C 3 , k m = C m /G. 
The elastic modulus is plotted versus temperature T in Figure 14, and the constants, k m , 
are depicted in Figures 15 and 16. The experimental data are approximated by the linear 
equations 

G = Go — G\T, k m = k m o + k m iT, (14) 

where Gj and k m j (j = 0, 1) are found by the least-squares technique. The quantities 
(z) = 4Z and r are plotted versus temperature in Figure 17 together with their mean 
values (z) = 1.35 and r = 0.14 s _1 . 



6 Discussion 

Figure 14 shows that G is independent of T in the sub-T g region, and it decreases with 
T above T g . This conclusion is in accord with the conventional assertion that the elastic 
modulus falls down when the temperature reaches the glass transition point. A decrease 
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in G with temperature is quite noticeable (the modulus is reduced by 32 %), but it is 
substantially less than that in polymeric glasses (where elastic moduli drop by several 
orders of magnitude as temperature passes through the glass transition point). 

Figures 15 and 16 reveal that the concentration of "viscous" CRR increases with 
temperature both below and above the glass transition point, whereas the fractions of 
"elastic" and viscoelastic relaxing units remain practically constant below T g . The latter 
means that only a portion of passive CRRs with negligible potential energies become 
activated in the sub-T g region. Above the glass transition point, k\ and k 2 substantially 
increase with temperature, while k^ pronounsly decreases. This reflects transformation of 
passive relaxing units into active "viscous" and "viscoelastic" CRRs induced by thermal 
fluctuations. 

Figure 17 demonstrates that the attept rate, r , and the average activation energy 
of CRRs, (z), are independent of temperature. The fact that the energy landscape of a 
glass is not affected by temperature appears to be rather natural jlOfl . To explain why T 
remains temperature-independent, we refer to the Narayanaswamy theory JX1|J . According 
to it, rearrangement is CRRs is governed by some effective temperature, T e g, 

1 a 1 — a ,„ 

^ = r + — (15 > 

where Tf is the Active temperature (that characterizes changes in the internal structure of 
a glass), and a G [0, 1] is a material constant. With an increase in temperature, the first 
term on the right-hand side of Eq. (15) decreases. This decrease is, however, compensated 
by an increase in the other term: given a duration of a dynamic test, structural recovery 
in a glass occurs more intensively at a higher temperature, which implies that a drop 
in the fictive temperature (that approaches the current temperature for an equilibrated 
glass) is more pronounced compared to that at a low temperature. 

It is worth to compare Eq. (8) with other micro-mechanical concepts in viscoelas- 
ticity of glasses. The ADWP (asymmetric double-well potential) model |T2j does not 



account for the presence of "viscous" CRRs that play the key role in the time-dependent 
response (Figure 15). The SGR (soft glassy rheology) model pi presumes strain-induced 



transformations of the energy landscape that are not observed in experiments (Figure 17). 



7 Concluding remarks 

Constitutive equations have been developed for the viscoelastic behavior of a silicate glass 
at isothermal uniaxial loading with small strains. A glass is treated as an ensemble of 
independent units, where relaxation of stresses is driven by cooperative rearrangement 
of clusters of atoms. The rearrangement process is thought of as a series of hops of 
CRRs (trapped at the bottom levels of their potential wells on the energy landscape) to 
higher energy levels. The ensemble is split into three groups of relaxing units: "elastic" 
CRRs, whose rate of rearrangement is substantially less that the minimal frequency, 
"viscous" CRRs whose rate of rearrangement noticeably exceeds the maximal frequency, 
and "viscoelastic" CRRs, whose rate of rearrangement is comparable with frequencies of 
dynamic tests. Explicit expressions have been derived for the storage and loss moduli. 
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These formulas are determined by 5 adjustable parameters that are found by matching 
experimental data for a multicomponent silicate glass at temperatures in the vicinity 
of the glass transition point. It is demonstrated that the attempt rate for hops and 
the distribution of potential wells with various depths are independent of temperature, 
whereas the instantaneous elastic modulus and the concentrations of "elastic" , "viscous" 
and "viscoelastic" relaxing units are strongly affected by temperature. With an increase 
in T, the fraction of "elastic" units decreases, and the concentrations of "viscous" and 
"viscoelastic" units increase, which is ascribed to thermally-induced activation of passive 
CRRs. 
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Figure 1: The loss tangent r = log tan 5 versus temperature T °C. Circles: treatment of 
observations. Solid line: approximation of the experimental data by a piecewise linear 
function. The vertical line indicates the glass transition temperature 




Figure 2: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency u rad/s. Symbols: experimental data (T = 500.0 °C). 
Solid lines: numerical simulation 
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Figure 3: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency ui rad/s. Symbols: experimental data (T = 526.5 °C). 
Solid lines: numerical simulation 




Figure 4: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency u rad/s. Symbols: experimental data (T = 552.0 °C). 
Solid lines: numerical simulation 
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Figure 5: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency ui rad/s. Symbols: experimental data (T = 577.0 °C). 
Solid lines: numerical simulation 
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Figure 6: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency u rad/s. Symbols: experimental data (T = 601.5 °C). 
Solid lines: numerical simulation 
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Figure 7: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency ui rad/s. Symbols: experimental data (T = 627.5 °C). 
Solid lines: numerical simulation 




Figure 8: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency u rad/s. Symbols: experimental data (T = 653.0 °C). 
Solid lines: numerical simulation 
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Figure 9: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency u rad/s. Symbols: experimental data (T = 678.0 °C). 
Solid lines: numerical simulation 




Figure 10: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency u rad/s. Symbols: experimental data (T = 706.5 °C). 
Solid lines: numerical simulation 
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Figure 11: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency ui rad/s. Symbols: experimental data (T = 727.0 °C). 
Solid lines: numerical simulation 




Figure 12: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency u rad/s. Symbols: experimental data (T = 749.5 °C). 
Solid lines: numerical simulation 
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Figure 13: The storage modulus G' GPa (unfilled circles) and the loss modulus G" GPa 
(filled circles) versus frequency uj rad/s. Symbols: experimental data (T = 760.0 °C). 
Solid lines: numerical simulation 




Figure 14: The elastic modulus G GPa versus temperature T °C. Circles: treatment of 
observations. Solid lines: approximation of the experimental data by Eq. (14). Curve 1: 
G = 25.23, G 1 = 0.0; curve 2: G = 83.48, d = 8.63 • 10" 2 
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Figure 15: The dimensionless parameter k\ versus temperature T °C. Circles: treatment 
of observations. Solid lines: approximation of the experimental data by Eq. (14). Curve 
la: k 10 = -2.09 • 10~ 2 , k u = 5.31 • 10" 4 ; curve lb: k 10 = -2.43, k u = 3.48 • 10~ 3 




Figure 16: The dimensionless parameters k 2 (unfilled circles) and k 3 (filled circles) versus 
temperature T °C. Symbols: treatment of observations. Solid lines: approximation of 
the experimental data by Eq. (14). Curve 2a: k 2 o = 3.80 • 10~ 2 , k 2 \ = 0.0; curve 2b: 
k 20 = -5.34, k 21 = 8.02 • 10~ 3 ; curve 3a: k 30 = 0.95, k 31 = 0.0; curve 3b: k 30 = 8.10, 
k m = -1.06 • 10- 2 



16 




500 



800 



Figure 17: The average activation energy per CRR (z) (unfilled circles) and the attempt 
rate T s' 1 (filled circles) versus temperature T °C. Symbols: treatment of observations. 
Curves 1 and 2: the mean values 
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